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The Fermi Large Area Telescope (LAT) is a powerful pulsar detector, as demonstrated by the over one hundred 
objects in its second catalog of pulsars. Pass 8 is a new reconstruction and event selection strategy developed 
by the Fermi-LAT collaboration. Due to the increased acceptance at low energy, Pass 8 improves the pulsation 
detection sensitivity. Ten new pulsars rise above the 5 sigma threshold and are presented in this work, as well 
as one previously seen with the former Pass 7 reconstruction. 

More than 60% of the known pulsars with spin-down power ( E ) greater than 10 36 erg/s show pulsations in 
gamma-rays, as seen with the Fermi Large Area Telescope. Many non-detections of these energetic pulsars 
are thought to be a consequence of a high background level, or a large distance leading to a flux below the 
sensitivity limit of the instrument. The gamma-ray beams of the others probably miss the Earth. The new 
Pass 8 data now allows the detection of gamma ray pulsations from three of these high spin-down pulsars, PSRs 
J1828—1101, J1831—0952 and J1837—0604, as well as three others with E > 10 35 erg/s. We report on their 
properties and we discuss the reasons for their detection with Pass 8. 


1. Introduction 

Since its launch in June 2008, the Fermi satellite 
has accumulated thousands of hours of observations 
of the sky. Events recorded by the LAT have very ac¬ 
curate (< 1 qs) timestamps derived from GPS clocks 
onboard the satellite. This precise timing associated 
with the well-known position of the spacecraft allowed 
the detection of 132 pulsars listed in the second pul¬ 
sar catalog (2PC) of gamma-ray pulsars 1.]. These 
objects can be classified in different categories, allow¬ 
ing population studies, as shown in Section [2] 

Pass 8 is a new reconstruction and event selection 
strategy developed by the Fermi -LAT collaboration. 
It allows better acceptance and sensitivity at low en¬ 
ergies compared to Pass 7 data (P7REP), as described 
in Section [3] Most pulsars have spectra that cut off 
around a few GeV and therefore have most of their 
flux at these low energies. As a consequence, we ex¬ 
pect more pulsar detections in the future thanks to 
the new Pass 8 reconstruction. Of the eleven new 
pulsars presented in Section [3] (see Figs. [L] [2] and 
Table HI), ten were not seen with Pass 7. These new 
detections bring the total number of gamma-ray pul¬ 
sars known at present to 163 (see Fig. [4j. If we also 
count radio millisecond pulsars discovered at the posi¬ 
tion of unidentified gamma-ray sources but for which 
an ephemeris reliable enough to phase-fold the LAT 
data is not yet available, over 200 gamma-ray pulsars 
are now known. 

The fraction of pulsars detected in gamma-rays in¬ 
creases with spin-down power ( E ), as shown in Fig. 
[3] However, different factors such as distance or low 
signal-to-background ratio can complicate their detec¬ 
tion. This is the case for the six energetic pulsars 


(.E > 10 35 erg s _1 ) presented here, which are located 
in the central regions of our Galaxy. We focus on these 
high spin-down power pulsars in Section |4j 


2. General properties of gamma-ray 
pulsars 

Gamma-ray pulsars can be divided in two main cat¬ 
egories. “Young” ones are produced after a supernova 
explosion of a massive star and “recycled” ones are 
old pulsars re-accelerated after the accretion of mat¬ 
ter from a binary companion. The latter objects ro¬ 
tate much faster than the young ones with a period 
of the order of 1 to 30 ms, and are called “millisecond 
pulsars” (MSPs). The two main categories are well 
separated in the P — P diagram represented in Fig. [4] 
The fraction of MSPs among the Fermi -LAT pulsars 
currently amounts to 43% while the fraction of MSPs 
among the total number of known pulsars is only of 
the order of 10%. 

When a pulsar is already known from radio or X-ray 
observations, the corresponding ephemeris is used to 
search for pulsations in the LAT data. This technique 
allowed the detection of more than a hundred of young 
and recycled pulsars. 

“Blind period searches” of unidentified LAT 
sources, in LAT data and with radiotelescopes, led 
to dozens of new pulsars. Deep radio follow-ups of 
the gamma-ray pulsars discovered in the LAT blind 
searches determined that nearly all are “radio-quiet” 
(Si 400 < 30 pJy, see Fig. 3 in [lj). 
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Figure 1: Gamma-ray (red) and radio (black) phase-folded light-curve of each new millisecond pulsar. The x-axis is the 
phase rotation of the pulsar and the y-axis corresponds to the weighted gamma-ray counts as obtained with the 
probability-weighting method described in the text. The radio data come from the Nangay radiotelescope [1] but PSR 
J0931—1902 was first detected with the Green Bank Telescope [§}. 


Geometry determines the radio and gamma-ray 
beam shapes. It depends on the angles of the mag¬ 
netic and rotation axes relative to each other and to 
the line-of-sight from the Earth. It also depends on 
how the “light cylinder” (radius r = cP/2ir , where 
an object in co-rotation with a neutron star turning 
with period P would reach the speed of light) crosses 
the not-quite-dipole shaped magnetic field. Therefore, 
beam shapes are observables that can constrain emis¬ 
sion models. Eerrm-LAT’s 40 radio-quiet pulsars are 
precious in this regard: nearly all other known pulsars 
are seen with radio telescopes (“radio-loud”), mean¬ 
ing that they all have geometries such that the ra¬ 
dio beam sweeps the Earth. The radio-quiet pulsars 
have the radio beams tilted elsewhere. The gannna- 
ray beams are very narrow in neutron star longitude 
(due to concentration of the gamma-radiating elec¬ 
trons and positrons along “caustically” focussed mag¬ 
netic held lines), but are very broad in latitude, being 
brightest near the neutron star equator, and fading to¬ 
wards the poles. Romani & Watters UM use these 
arguments to generate an “Atlas” of gamma-ray beam 
profiles, including tallies of the numbers of radio-loud 
versus quiet pulsars. The Eerrm-LAT pulsars also dif¬ 
fer from those often chosen for geometry studies by 
their large spin-down powers indicative of strong mag¬ 
netic fields and relatively short periods, resulting in 
different magnetospheric configurations. 


3. New detections with Pass 8 

In the beginning of the mission, event reconstruc¬ 
tion was based on pre-launch instrument simulations. 
This reconstruction was close to reality but after 
analysing the first data, the Eermi-LAT collaboration 
realized that due to residual signals induced by back¬ 
ground cosmic-rays the selection efficiency was not op¬ 
timal, in particular at lower energies. Therefore the 
simulations were improved in order to take into ac¬ 
count this effect and to better characterize the instru¬ 
ment performance. A completely new reconstruction 
was developed to limit the loss of data due to parasite 
signals. The event selection is now optimized and the 
systematic errors are significantly reduced. Together 
with this new reconstruction called “Pass 8” Q, the 
collaboration produced corresponding diffuse models 
to describe the extragalactic diffuse emission as well 
as gamma-ray emission due to cosmic-ray interaction 
with the Earth’s atmosphere. The Galactic diffuse 
model from the previous reconstruction P7REP was 
scaled to take into account the enhancement of the 
emission expected from Pass 8 acceptance improve¬ 
ment. This new reconstruction shows a gain in ef¬ 
fective area of 30% above 1 GeV and a factor 2 at 
100 MeV compared to the previous reconstruction 
P7REP. The angular resolution is also improved, lead¬ 
ing to a sensitivity gain of 40% for point-like sources. 

So far ten new pulsars have indeed been detected 
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Figure 2: Gamma-ray (red) and radio (black) phase-folded light-curve of the newly detected young pulsars. The radio 
data come from the Parkes telescope [fa] for PSRs J1224—6407, J1739—3023, J1828—1101, J1831—0952 and 
J1837—0604 and from Jodrell Bank Observatory 0 for PSRs J1856+0113 and J1857+0143. 


exclusively with Pass 8. We also analysed one more 
pulsar, PSR J1856+0113, that was already seen with 
P7REP data (just above the 5er threshold) but for 
which Pass 8 improved the pulsation detection to over 
8cr. Hence we present here 11 new pulsars analysed 
with the Pass 8. Their gamma-ray and radio light- 
curves are presented in Fig. [T] (MSPs) and Fig. [2] 
(young) and their properties are listed in Table |I] 

To detect pulsations, regions of 3° around the ra¬ 
dio position of the pulsars were selected. The data 
were phase-folded using radio eplremerides from differ¬ 
ent radiotelescopes: Nangay [dj, Parkes |[6j and Jodrell 
Bank |7). A probability-based event selection was then 
applied, using the shape of the point-spread function 
to estimate the events’ probability of originating from 
the pulsar. The lowest H-test value obtained for this 
sample of pulsars is ~40 for PSR J1828—1101. 

The spectral analysis was performed using the 
Fermi ScienceTools and after selecting a region of 


15° around each pulsar. We restricted to the energy 
range between 100 MeV and 300 GeV and we used the 
Source class events. Sources from the third Fermi- 
LAT catalog (3FGL) [H were included in the model 
of the regions and a point-like source with a power-law 
spectral shape was added at the center of each ROI, 
corresponding to the position of each pulsar. PSRs 
J0742+6620, J0931-1902 and J1837-0604 were coin¬ 
cident with unidentified 3FGL sources, which we as¬ 
sumed to stem from the pulsars, therefore the position 
of the corresponding source was shifted to the radio 
position of the pulsar. When performing the fit of the 
region with gtlike, the spectral parameters of sources 
within a radius of 5° were let free as well as the ones 
within 10° and with a TS value above 1000 (taken ar¬ 
bitrarily) in the 3FGL. No phase selection was applied 
to the data for this analysis. 

Although all pulsars presented here are significantly 
detected with their pulsations in gamma-rays, their 
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spectral analysis was not successful in many cases. 
For all but two of the pulsars (J1730—2304 and 
J1857+0143), the light-curves show very narrow peaks 
(see Figs. |T] and [2]) easily detected by phase-folding. 
However when considering the full phase-band for the 
spectral analysis, the faint signal fades behind the 
background level and the source is not detected, as 
can be seen by the TS value column in Table |T) An 
analysis selecting only the on-phase intervals for each 
pulsar will be made for future publication. 

Among the 11 new objects, 4 are MSPs and two of 
them have a period of ~ 8 ms. These two new detec¬ 
tions start to fill the bridge between MSPs and young 
populations in Fig. [I] All the new MSPs but one 
(J1730—2304) are located far from dense background 
regions (with a latitude |6| > 20°) where we can de¬ 
tect very faint objects such as J0931 —1902 which is 
the pulsar with the lowest energy flux measured with 
the LAT at present (see Table Q] as well as 0 , 0 )- 

4. Focus on high spin-down power 
pulsars 

More than 60% of pulsars with E > 10 36 erg/s 
are detected in gamma-rays, and more than 50% for 
E > 10 35 erg/s, as can be seen in Fig. [3] Six of the 
new gamma-ray pulsars presented here belong to the 
energetic pulsar category, E > 10 35 erg/s. 

The pulsar timing campaign for Fermi |9] focused 
on pulsars with E > 10 34 erg/s. In Figure [3] (top) 
nearly 90% of the young high E pulsars have indeed 
been monitored. The lower rate for MSPs is due to 
those discovered after the campaign list was estab¬ 
lished. The fraction is lower for low E pulsars, but 
they are abundant, resulting in a large sample never¬ 
theless. The choice of which low E pulsars are mon¬ 
itored by radio telescopes could conceivably lead to 
bias in the fraction of gamma-detected pulsars shown 
here. Pass 8 makes the LAT more efficient at finding 
pulsars and thus reduces the biases in determining 
these fractions. 

A very large fraction of high E MSPs are gamma- 
ray pulsars: the small light-cylinder leads to very 
broad beams. See M for a discussion of those not 
seen in gamma-rays. The gamma-ray deathline is at 
lower E for MSPs compared to slower pulsars. For 
the young pulsars, the fraction increases with E , pass¬ 
ing the 50% mark around 10 35 erg/s. Luminosity in¬ 
creases with E, but sensitivity and background levels 
account for only part of the undetected pulsars. Beam 
geometry is surely the dominant factor: the “Atlas” 
of gamma-ray profiles provided by 0 , and similar 
studies since then, show that models do indeed predict 
roughly that fraction of pulsars where the gamma-ray 
beam either misses the Earth, or is so broad as to 
give unpulsed detection. Hou et al. 0] discuss this 



logA (erg s 



I og ( 'J E j't-Kd 2 ) ( \J erg s 1 kpc 2 ) 


Figure 3: For each decade in spin-down power (top) or in 
heuristic gamma-ray flux (bottom, see text), the solid 
green line shows the fraction of known “young” pulsars 
(Po > 30 ms) for which we have rotation ephemerides, 
and the dashed red line shows the same for MSPs (Po < 
30 ms). We consider only pulsars outside of globular 
clusters. We have gamma-ray phase-folded all pulsars for 
which we have ephemerides. The green dots (red 
triangles) show the fraction of these gamma-ray 
phase-folded young (millisecond) pulsars for which the 
LAT sees pulsations. The error bars are the 68% 
confidence level statistical uncertainties on the fraction, 
using the Bayesian calculation of El¬ 


is more detail. 

The “heuristic” gamma-ray flux, \/Tj/d 2 (see 2PC 
equation 18 ) uses the idea that gamma-ray luminos¬ 
ity L 1 scales with the open field line voltage oc x/lt’, 
loosely born out by the correlation between L 7 and 
E seen in 2PC Figure 9. The solid line in Figure [3] 
bottom, differs from that in the top frame because of 
the radio-quiet pulsars for which we have no distance 
estimate. Figure [3] shows (as does 2PC Figure 15) 
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Table I : Temporal and spectral properties of the new pulsars detected with Pass 8. The first four lines are the new 
MSPs. The period and spin-down power are taken from 0. The distance in Column 4 is estimated from the 
dispersion measurement using the NE2001 model 0, and the uncertainties come from re-running NE2001 for ±20% of 
the DM. Column 5 corresponds to the test statistic value obtained after the spectral analysis described in the text (*), 
or after a study with the pointlike tool(+). Gioo is the integrated energy flux between 0.1 and 300 GeV assuming a 
power-law spectrum with the corresponding index value —T. Columns 8 and 9 give the total gamma-ray luminosity in 
the 0.1 to 300 GeV energy band and the gamma-ray conversion efficiency p = L^/E. The first uncertainty in L 7 and p 
comes from the statistical uncertainties in the spectral fit while the second is due to the distance uncertainty. When 
the pulsar is too weak or lying in background-contaminated regions, we could not derive a spectrum, therefore we do 
not give any spectral information. 


PSR Name 

Period 

E 

Distance 

TS 

r 

Gioo 

l 7 

Efficiency 


( ms ) 

(erg s -1 ) 

(kpc) 



(10” n erg cm" 2 s 1 ) 

(10 33 erg s 1 ) 

(%) 

J0742T6620 

2.89 

2.0xl0 34 

0.68 ±0.10 

143* 

2.4±0.1 

0.49 ± 0.06 

0.27 ±0.03 ±0.07 

1.3 ±0.2 ±0.4 

J0931—1902 

4.64 

1.4x 10 33 

1.88 ±0.51 

23* 

2.0±0.2 

0.22 ±0.06 

0.9 ± 0.2^ 4 

65 ± 18t 46 

J1455—3330 

7.99 

1.9x 10 33 

0.53±0.07 

0.5+ 

- 

- 

- 

- 

J1730—2304 

8.12 

1.5xl0 33 

0.53±0.05 

65* 

2.6±0.05 

1.16 ±0.09 

0.39 ±0.03 ±0.07 

26 ± 2 ± 5 

J1224—6407 

216.50 

1.9x 10 34 

3.15 ± 0.41 

0.04+ 

- 

- 

- 

- 

J1739—3023 

114.37 

3.0x10 s5 

2.91 ±0.38 

29* 

2.33±0.005 

1.61 ±0.02 

16.2 ±0.2 ±4 

5.27 ± 0.07 ± 1 

J1828—1101 

72.05 

1.6x 10 36 

6.63±1.05 

73* 

2.5±0.1 

2.7 ±0.4 

140 ± 20 ± 40 

9 ± 1 ± 3 

J1831—0952 

67.27 

l.lxlO 36 

4.05±0.37 

2.3+ 

- 

- 

- 

- 

J1837—0604 

96.29 

2.0xl0 36 

6.41±0.67 

439* 

2.50±0.05 

7.6 ±0.5 

370 ± 30 ± 80 

19 ±1 ±4 

J1856T0113 

267.44 

4.3xl0 35 

3.07±0.32 

0.4+ 

- 

- 

- 

- 

J1857T0143 

139.76 

4.5x 10 35 

5.75±0.44 

0.0+ 

- 

- 

- 

- 


that for \/H/d 2 below 10 15 -\/erg s -1 kpc 2 the num¬ 
ber of detected pulsars falls to zero. This is due both 
to the LAT’s flux sensitivity and to the E deathline 
near 10 33 erg s _1 seen in both Figure[4]and the top of 
Fig. H 

We note that all of the newly detected energetic 
pulsars lie in very crowded regions close to the galac¬ 
tic center and are therefore subject to a high diffuse 
emission level. Three of them are also quite far away 
with a distance estimate above 5 kpc. Pass 8 demon¬ 
strates here its ability to detect faint sources above 
the background, with its sensitivity gain of 40% for 
point-like sources. 

Finally, an important factor for pulsation detection 
is also the quality and completeness of the radio tim¬ 
ing which can be achieved thanks to the precious col¬ 
laboration of radiotelescopes teams. 


Conclusion 

We presented 4 new MSPs and 7 new young pulsars 
detected in gamma-rays with the Fermi- LAT. Two 
of the MSPs are among the faintest gamma-ray pul¬ 
sars detected at present with an energy flux below 
5 x 10~ 12 erg cm' 2 s _1 . Six of the new young pulsars 
are very energetic objects, with E > 10 35 erg s _1 , and 
located in the central regions of the Galaxy. Among 
the undetected energetic pulsars, different limitations 
prevent the detection such as a large distance induc¬ 
ing a flux below the sensitivity limit of the instru¬ 


ment; a high background level leading to a low signal- 
to-noise ratio; or intrinsic pulsar features (low cutoff, 
wide peaks, beam sampling...). There is no hope to 
detect pulsars whose gamma-ray beam do not sweep 
Earth, but the other limitations can be overtaken with 
the increased acceptance of Pass 8 data which will 
certainly help detecting more new objects, as it was 
demonstrated in this work. 
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Figure 4: P -P diagram with red triangles representing the gamma-ray MSPs, blue squares are the radio-quiet 
gamma-ray young pulsars and green dots are radio-loud gamma-ray young ones. We emphasized the newly detected 
objects. Black dots are timed but undetected pulsars in gamma-rays and grey ones are pulsars we have not 
phase-folded. 


References 

[1] Abdo, A. A. et al., “The Second Fermi Large 
Area Telescope Catalog of Gamma-Ray Pulsars”, 
ApJS, 208, 17 (2013) 

[2] Romani, R. W. and Watters, K. P, “Constraining 
Pulsar Magnetosphere Geometry with gamma- 
ray Light Curves”, ApJ, 714, 810 (2010) 

[3] Atwood, W. et al., “Pass 8: Toward the 
Full Realization of the Fermi-LAT Scientific 
Potential”.arXiv:1303.3514 (2013) 

[4] Cognard, I. et al., “Discovery of Two Millisecond 
Pulsars in Fermi Sources with the Nangay Radio 
Telescope”, ApJ, 732, 47 (2011) 

[5] Lynch et al., in preparation 

[6] Weltevrede, P. et al., “Pulsar Timing with the 
Parkes Radio Telescope for the Fermi Mission”, 
PASA, 27, 64 (2010) 

[7] Hobbs, G. et al., “Long-term timing observations 
of 374 pulsars”, MNRAS, 353, 1311 (2004) 

[8] Acero, F. et al., 11 Fermi Large Area Telescope 
Third Source Catalog”, arXiv:1501.02003 (2015) 

[9] Smith, D. A. et al., “Pulsar timing for the 


Fermi gamma-ray space telescope”, AAP, 492, 
923 (2008) 

[10] Guillemot, L. and Tauris, T. M., “On the non¬ 
detection of gamma-rays from energetic millisec¬ 
ond pulsars - dependence on viewing geometry”, 
MNRAS, 439, 2033 (2014) 

[11] Watters, K. P. et al., “An Atlas for Interpreting 
gamma-ray Pulsar Light Curves”, ApJ, 695, 1289 
(2009) 

[12] Hou X., et al., “Six faint gamma-ray pulsars seen 
with the Fermi Large Area Telescope. Towards a 
sample blending into the background”, AAP, 570, 
A44 (2014) 

[13] http://www.atnf.csiro.au/people/pulsar/psrcat/ 

[14] Cordes, J. M. and Lazio, T. J. W., “NE2001.I. 
A New Model for the Galactic Distribu¬ 
tion of Free Electrons and its Fluctuations”, 
astro-ph/0207156 (2002) 

[15] Paterno, M., “Calculating Efficiencies and Their 
Uncertainties” , FERMILAB-TM-2286-CD, May 
22, 2004. 


eConf C141020.1 














